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Abstract 

One of several possibilities to construct a quantum theory of gravity is 
employing the Feynman path integral. This approach is plagued by some 
problems: the integration measure is not uniquely defined, the Einstein- 
Hilbert action unbounded, and perturbation theory nonrenormalizable. 
To make the path integral tractable one can approximate the continu- 
ous geometry of spacetime by a simplicial complex. The edge lengths of 
this lattice are considered as the dynamical degrees of freedom and Regge 
calculus is applied. In this work, numerical simulations using the Regge- 
Einstein action and a "compact" action show the occurence of a phase 
transition. The strength of this transition, separating a well-defined phase 
with finite expectation values from an ill-defined phase, is weaker for the 
compact action, which might be important for the continuum limit. To 
analyze the interaction mechanism of this formulation of quantum grav- 
ity, correlation functions between geometrical quantities like edge lengths, 
volume elements, and local curvatures have been computed. Our results 
for the two-point functions seem to prefer exchange particles with an ef- 
fective mass. To ease treatment of quantum gravity a new approach is 
proposed consisting in a transformation of the path integral to the parti- 
tion function of a spin system. This facilitates analytical and numerical 
calculations considerably. First results for the phase structure in two as 
in four dimensions are presented and indicate promising similarities to the 
original Regge theory. 
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1 Technical and Conceptual Problems of 
Quantum Gravity 



At present, with quantum chromodynamics for the strong interaction, the 
Glashow-Salam- Weinberg model for the electromagnetic and weak interactions, 
and Einstein's theory of general relativity for gravitation one may feel to have 
arrived at a closed system of physical laws. However, problems arise when trying 
to describe these theories in a unified framework. Although particle physics is 
successfully described by quantum field theories, the laws of gravity determining 
the large scale structures of the universe seem to be different. Based on the 
assumption that quantum theory is universally valid it would appear strange to 
have a drastically different framework for gravitation. Notwithstanding many 
decades of intense work, there is still no complete quantum theory of gravity 
with the properties thought to be essential for consistency, such as unitarity, 
renormalizability (or some related property, like finiteness or asymptotic safety) 
and Lorentz invariance in a local inertial frame. The difficulties are compounded 
by the total lack of any empirical data that is manifestly relevant to the problem 
and a serious obstacle for experiments in this field is the extreme smallness of 
the Planck length 



There are three especially important issues confronting every approach to quan- 
tum gravity [f, 2, 3]. 

1.1 Background structure 

The mathematical model of spacetime used in classical general relativity is a 
differentiable manifold equipped with a Lorentzian metric. Now one can ask for 
the underlying substructure of this picture. 

The bottom level is a set M whose elements are to be identified with spacetime 
points or events. This set is formless with its only general mathematical property 
being the cardinal number. In particular, there are no relations between the 
elements of M and no special way of labeling any such element. The next step is 
to impose a topology on M so that each point acquires a family of neighborhoods. 
It now becomes possible to talk about relationships between points, albeit in a 
rather nonphysical way. This defect is overcome by adding the key ingredient 
of all standard views of spacetime: the topology of M must be compatible with 
that of a differentiable manifold. A point can then be labeled uniquely in M 
(at least locally) by giving the values of n real numbers, with n the dimension 
of the manifold. Such a coordinate system also provides a more specific way of 
describing relationships between points of Af , though not intrinsically in so far 
as these depend on which coordinate systems are chosen to cover M. In the final 
step a Lorentzian metric g is placed on Af , thereby introducing the ideas of the 
length of a path joining two spacetime points, parallel transport with respect to 
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a Riemannian connection, causal relations between pairs of points etc. 

However, since one wishes to assert that some sort of quantum spacetime 
structure is meaningful, the key question for any particular approach to quantum 
gravity is how much of the hierarchy described above must be kept fixed. For 
example, in most of the former approaches to quantizing the gravitational field, 
the set of spacetime points, topology and differential structure are all fixed, and 
only the Lorentzian metric g is subject to quantum fiuctuations. If the Lorentzian 
metric g becomes quantized then the light cone associated with any spacetime 
point is no longer fixed and it is not meaningful to impose microcausal relations. 
This destroys one of the bedrocks of conventional quantum field theory and is 
probably the greatest reason why spacetime approaches to quantum gravity have 
not got as far as might have been hoped. 

1.2 Spacetime difFeomorphism group 

The group Diff(Af) of spacetime diffeomorphisms plays a key role in the clas- 
sical theory of general relativity and so the question of its status in quantum 
gravity is of considerable interest. The action of Diff(Af) on M affects the space 
of spacetime fields, and the only thing that has immediate physical mean- 
ing is the quotient space J-'/Diff(Af) of orbits, i.e. two field configurations are 
regarded as physically equivalent if they are connected by a Diff(Af) transforma- 
tion. Technically, this is analogous to the situation in electromagnetism where 
a vector potential is equivalent to -|- 9^/ for all functions /. However, 
there is an important difference between electrodynamics and general relativity. 
Electromagnetic gauge transformations occur at a fixed spacetime point x, and 
the physical configurations can be identified with the values of the electromag- 
netic field Ff^i,{x), which depends locally on points of M. On the other hand, 
Diff(Af ) maps one spacetime point into another, and therefore the obvious way of 
constructing a Diff(Af )-invariant object is to take some scalar function of space- 
time fields and integrate it over the whole of Af , which gives something that is 
very nonlocal. The idea that physical observables are naturally nonlocal is an 
important ingredient in many approaches to quantum gravity. 



One of the major issues in quantum gravity is the so-called problem of time. This 
arises from the very different roles played by the concept of time in quantum 
theory and in general relativity. In standard quantum theory time is not a 
physical observable in the normal sense since it is not represented by an operator. 
Rather, it is treated as a parameter which, as in classical physics, is used to mark 
the evolution of the system. In particular, it provides the parameter t in the 
time-dependent Schrodinger equation 



1.3 Time 



= Hi^ ■ 

dt ^ 
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In addition, the scalar product of states and complete sets of commuting observ- 
ables are defined at a single time, measurements are carried out at a single time, 
and states specify the probability of the results of such measurements. 

In a Diff(Af )-invariant theory like classical general relativity the role of time is 
very different. If M is equipped with a Lorentzian metric </, and if its topology is 
appropriate, it can be foliated in many ways as a I-parameter family of spacelike 
surfaces, and each such parameter might be regarded as a possible definition of 
time. Such a definition of time is rather unphysical since it provides no hint as 
to how it might be measured or registered. Moreover, the possibility of defining 
time in this way is closely linked to a fixed choice of the metric g. It becomes 
untenable if g is subject to some type of quantum fiuctuation. 

One might try to construct a theory of quantum gravity by using functional 
integrals analogous to those employed in normal quantum field theory to produce 
the vacuum expectation values of the time-ordered product of fields. However, 
it is difficult to see what a time-ordered product could mean in the absence of 
any background metric to provide a preferred notion of timelike and spacelike. 

Attempts to construct a quantum theory of gravity can be divided into three 
broad categories briefiy described in the following Sections. 

1.4 Quantizing general relativity 

The idea is to start with the classical theory of general relativity represented by 
the Einstein-Hilbert action 

[G labels Newton's constant, g the determinant of the metric tensor, and R 
the curvature scalar) and then to apply some type of quantization algorithm. 
This is intended to be analogous to the way in which the classical theory of 
an electron bound by the Coulomb potential is quantized by replacing certain 
classical observables with self-adjoint operators on a Hilbert space. Of course, 
this is essentially also the approach used in developing important elementary- 
particle physics ideas like the Glashow-Salam- Weinberg electro-weak theory and 
the quantum chromodynamics description of the strong nuclear force. 

In early particle-physics based approaches to quantum gravity the starting 
point is to fix the background topology and differential structure of spacetime 
M to be that of Minkowski space, and then to write the Lorentzian metric g on 
M as 

gab{x) = rjab + hab{x) , (4) 

where h measures the departure of g from fiat spacetime rj. The background met- 
ric J] provides a fixed causal structure with the usual family of Lorentzian inertial 
frames. Thus, at this level, there is no problem of time. The causal structure 
also allows a notion of microcausality, thereby permitting a conventional type of 
relativistic quantum field theory to be applied to the field hah- The quanta of 
this field are massless spin-2 particles, the gravitons. A typical task would then 
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be to compute perturbative scattering-matrix elements for these gravitons. The 
action of Diff(Af) is usually studied infinitesimally and is reflected in the quan- 
tum theory via a set of Ward identities that must be satisfied by the n-point 
functions of the theory. Unfortunately, ultraviolet divergences are sufficiently 
violent to render the theory perturbatively nonrenormalizable [4]. 

One of many responses was to enlarge the classical theory of general relativity 
with carefully chosen matter fields with the hope that the ultraviolet divergences 
would cancel, leaving a theory that is perturbatively well-behaved. The concept 
is based on the extension of the group of global Poincare transformations with 
anticommuting spinor generators Q. Hence, a supersymmetry between bosons 
and fermions is achieved, which allows to transform bosonic to fermionic states 
and vice versa 

Q\B) ~ \F) , Q\F) ~ \B) . (5) 

Thus, each particle obtaines a corresponding SUSY partner. Indeed, diver- 
gent bosonic (fermionic) loop contributions can compensate those of fermions 
(bosons). A generalization to local supersymmetric transformations leads to a 
gauge theory of supersymmetry, called supergravity [5] . The Lagrangian of sim- 
ple supergravity reads 

C = --^eR{e,co) - ie'='™-Vfc757/^™V'n , (6) 
IdttG z 

with e the determinant of the vielbein, io the spin connection, and 7i the Dirac 
matrices. The gravitino field is the SUSY partner of the graviton and is 
associated to the generator Q. Early expectations were following successful low- 
order results but it is now generally accepted that if higher-loop calculations 
could be performed intractable divergences would appear once more. 

The fact that a theory is perturbatively nonrenormalizable does not neces- 
sarily mean that it has an intrinsic problem or that it is bad in any way. It 
merely says that perturbation theory does not apply to the problem in ques- 
tion. It is therefore not surprising that a great deal of developments focus on 
nonperturbative quantization in the context of the canonical theory. 

The first step in canonically quantizing general relativity consists in foliating 
spacetime into a family of space-like 3-dimensional hypersurfaces, i.e. one decom- 
poses spacetime into space and time. The metric on these hypersurfaces, /^^^(x), 
will play the role of the canonical variable. All quantities are then decomposed 
into variables which exist on such a hypersurface and variables which point in 
the fourth, timelike dimension. It then turns out that the canonical momentum, 
n"^(x), is given by the extrinsic curvature of a 3-dimensional hypersurface, i.e. 
the quantity which describes the embedding of space into spacetime. As a con- 
sequence of invariance under arbitrary coordinate transformations one finds that 
there exist four constraints, of which the so-called Hamiltonian constraint 7i is 
directly connected to the invariance of the theory under reparametrizations of 
time. Its explicit form is 

n = ^^Gabcdii^'m' - T^VhR + nm = o, (7) 
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where R is the curvature scalar on 3-space and Tim is the Hamiltonian density 
for nongravitational fields. The coefficients Gabcd depend explicitly on the metric 
and play the role of a metric in the space of all metrics. Quantization now 
proceeds, by elevating the metric and its momentum to the status of operators 
and imposing the commutation relations 

[Kb{^),W'{y)] = zn6:6,'6{^-y) . (8) 

One specific realization of (8) is provided by the substitution 

^ -JJ- . (9) 

The classical constraint (7) is then formally implemented in the quantum theory 
by inserting (9) into (7) and applying it on wave functionals ^ which depend on 
the 3-metric and on nongravitational fields denoted by ^, i.e. 

nnKbi^), <^(x)] = (-^^^^Gabcd^^-r- - T^^R + nm]^ = o. (to) 

\ ohabhcd IbirG J 

This is the well-known Wheeler-DeWitt equation [6], the gravitational analog of 
the Schrodinger equation. The quantization of the remaining constraints leads 
to the condition that this wave functional does not change under coordinate 
transformations of the 3-metric, but is a function of the geometry only. The 
configuration space is thus the space of all 3-geometries and is called superspace. 
Several techniques have been suggested for finding solutions to the Wheeler- 
DeWitt equation. An interesting development has been the discovery of appro- 
priate variables which enables to find exact solutions of (7) in the absence of 
matter. This became possible since the complicated potential term of (7) disap- 
pears when it is rewritten in terms of these Ashtekar variables [7], which have a 
strong similarity to Yang-Mills loop variables. The solutions can be classified in 
terms of loops and knots and exhibit an interesting structure of space, of which 
the most important may be the existence of a minimal length [8]. 

Another nonperturbative scheme is provided by the functional integral ap- 
proach. Here only a brief overview is given, but it will be outlined in more 
detail in the next Chapter. In ordinary quantum mechanics a solution to the 
time-dependent Schrodinger equation (2) can be generated in the form 



(11) 

where the integral is over all paths that end at the point x at time t. The solution 
thus obtained depends on the value x[to) of the path at the initial time to. 

The gravitational analog of (II) plays an important role in the Hartle- 
Hawking approach to quantum gravity [9]. There, the wave functional is ex- 
pressed formal path integral 

^h] = I D[g]e-'-^^^ , (12) 

rJM 
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where the sum is over Euchdean geometries g on the 4-manifold M . The bound- 
ary dM = {Ei, Ej} consists of 3-dimensional hypersurfaces with induced metrics 
h. IeIq] is the Euchdean form of (3) and contains an additional boundary term. 
Note, that there is no explicit notion of time in (12) and the proper time between 
the surfaces depends on the 4-geometries in the sum. 

When applied to cosmology, ^ is the wave function of the universe. The main 
motivation for this claim comes from the nonseparability of quantum theory, i.e. 
from the fact that one cannot in general assign a wave function to a given system 
because it is not isolated but coupled to its natural environment, which is again 
coupled to another environment, and so forth. The extrapolation of this quantum 
entanglement leads inevitably to the concept of a wave function for the universe. 

It can be shown that (12) satisfies the Wheeler-DeWitt equation and again 
the solution depends on the boundaries. In the Hartle-Hawking no boundary 
ansatz it is assumed that one has to perform a sum over compact manifolds with 
one boundary only, which is given by the considered universe. The lack of a 
second boundary (at a small size of the universe) saves one from the need to find 
appropriate boundary conditions there: the universe is created ex nihilio. 

However, the quantum-gravity path-integral is also plagued with a lot of 
problems concerning the integration over 4-geometries, the Euclidean gravita- 
tional action and the measure. In particular, the integration should include a 
summation over manifolds of different topologies, but topologies are not classifi- 
able in 4 dimensions. Due to conformal fiuctuations, the gravitational action is 
not bounded from below, which results at a first view in a divergent path inte- 
gral. A unique definition of the measure does not exist because different physical 
motivations lead to different measures. 

Nevertheless, approximate calculations of functional integrals have been ex- 
plored to some extent. These include minisuperspace models obtained by freezing 
all but a finite number of modes of the 4-geometry, e.g. lattice gravity approaches 
like Regge calculus [10]. 

The idea of approximating spacetime by a simplicial complex has been of 
interest in both classical and quantum gravity for a long time. The Regge calculus 
provides a direct route in approximating continuous 4-geometries by simplicial 
lattices [11]. Expectation values can be written as 

{O) = Y.v{K) j dij\l]0{l)e-'-^\ (13) 

where ^{K) is the measure weight-factor of the simplicial complex K with A^o 
vertices, //[/] is the measure and 0{l) is an observable written as a function of 
the edge lengths / of the simplices. The Regge analog Ir of the Einstein-Hilbert 
action (3) is usually given as a function of the area At of the 2-simplices t and 
of the deficit angles 8t of the 4-simplices containing these faces 

Hl) = ^.^AA. (14) 

As the notation suggests it can be rewritten as a complicated function of the edge 
lengths alone. It is hoped that the classical continuum limit will be achieved as 
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A^o — ^ and / — )• 0. The other possibihty would be a finite length scenario with 
the attractive features of avoiding difficulties with locality at the Planck length 
and of providing a natural ultraviolet cut-off making renormalization of infinities 
unnecessary [12]. But still, many of the intriguing questions and problems arising 
in the conventional (continuum) functional integral approach are also present in 
the Regge calculus. 

1.5 "General-relativize" quantum theory 

The main idea is to begin with some prior concept of quantum theory and then 
to force it to be compatible with general relativity. The biggest programme of 
this type is due to Fredenhagen and Haag [f3]. It can be shown that the basic 
tenets of the general theory of relativity, namely general covariance and strict 
locality can be incorporated into quantum theory. Using the assumption that 
the theory has a scaling limit, each quantum state defines a reduced theory in 
the tangent space of each point. This reduced theory turns out to be invariant 
under translations and even under SLR{^). Then the macroscopic metric should 
evolve as a cooperative effect in finite size regions. However, many questions 
remain to be answered before a viable theory can be proposed. 

1.6 General relativity as a low-energy limit 

Another possibility is to remain with the idea of quantizing a classical system 
but start with something else than general relativity. The key step here is to 
find a system whose quantum theory is well-defined and which yields classical 
general relativity as a low-energy limit. A very sophisticated example for a 
scheme like this is the string theory, which abandons the concept of idealized 
point particles. The fundamental entities are I-dimensional objects (strings) 
instead of local quantum fields and the graviton occurs as just one of an infinite 
number of particles associated with the quantized string [14]. 
In the bosonic case one has the Polyakov action 

I[g. X] = J^J d'a./^)q^^ [a] d,X^ [a] d,X^ [a] g,^ [X] , (15) 
w 

where metric on the 2-dimensional manifold W (the world sheet) 

parametrized by cr = ((71,(72), and X : W ^ M are the bosonic string fields 
which map W into the spacetime manifold M (the target space) with back- 
ground metric (/^j,. The Regge slope a' is related to the string tension T = 
and is assumed to be of the order of the Planck length. 

A supersymmetric extension, the so-called heterotic string theory is formu- 
lated in a lO-dimensional world and obviously fails to reproduce the fact that we 
live in a 4-dimensional almost fiat Minkowski spacetime. Thus lO-dimensional 
superstring theory has to be compactified: six coordinates are curled up describ- 
ing a tiny compact space whose size is of the order of the Planck length, whereas 
the remaining physical four coordinates are kept uncompactified. Although it 
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was recently discovered how to construct string theories using conlormal field 
theory techniques without referring to any compactification scheme, there exists 
so far no compelling principle to determine the number of spacetime dimensions 
to be four. All dimensions below ten seem to be on an equal footing. 

Superstring theory has the great advantage over the covariant approaches 
that the individual terms in the appropriate perturbation expansion can be finite 
and, furthermore, the particle content of theories of this type could well be such 
as to relate the fundamental forces in a unified way. The low-energy limit of 
these theories is a form of supergravity because, when the string tension tends 
to infinity [a' — )• 0) the strings degenerate into points. Nevertheless, standard 
spacetime ideas in the sense of general relativity do not play a very significant 
role. This is refiected by the graviton being only one of an infinite number of 
particles in the theory. Similarly, the spacetime diffeomorphism group appears 
only as part of a much bigger symmetry structure. 
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2 Sum-over-Histories Approach 



In this work a nonperturbative treatment of quantum gravity using the path 
integral is persued. The essential point of this approach is Feynman's idea to 
represent the transition amplitude from an initial state with metric hi on a 
surface Ei to a final state with h2 on E2 as a sum over all field configurations [g] 
which take the given values on the surfaces Ei and E2 [3] 

{h2\hi) = f D[g]e^'^^^ , (16) 

M 

cf. Figure 1. The integration extends over all 4-geometries with the measure 
D[g]^ and I[g] denotes the gravitational action, which in general relativity is 
usually chosen as 

^=T^< [d'^V^{R-2A). (17) 
IdttG J 

M 

R is the curvature scalar, A the cosmological constant, and g the determinant of 
the metric. G is Newton's constant and units are such that c = h = 1. 

Strictly speaking the gravitational action contains second derivatives of the 
metric which have to be removed by integration by parts to give an action 
quadratic in first derivatives as is required by the path integral approach. There- 
fore, one gets a surface term 

/ d^xV±hK + C , (18) 

dM 

where K is the trace of the extrinsic curvature of the boundary dM with the 
induced metric h. The plus or minus sign is chosen according to whether the 




Figure 1: Manifold M with metric g and the boundary dM consisting of Ei and 
E2 on which the metrics hi and h2 are induced. 
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boundary is spacelike or timelike. This term turns out to be important to fulfill 
the composition law for the amplitude to go from the initial to the final state, 
to be obtained by summing over all states on some intermediate surface E3 with 
metric (see Figure 1) 

{h2\hi) = J2{h2\hs){hs\h^) . (19) 

hs 

Under variations of the metric the action (17) is stationary if the Einstein equa- 
tions in vacuum ^ 

Rmn--gmnR + Ag^^ = 0, (20) 

with Rmn the Ricci tensor, are satisfied. 

2.1 Gravitational action 

For real Lorentzian metrics g the action /[</, ^] will be real and so the integrand 
in the path integral will oscillate. In ordinary quantum field theory in fiat space- 
time one can circumvent this difficulty by performing a Wick rotation into the 
complex plane. It seems reasonable to apply similar ideas to the gravitational 
field. The analog of the Wick rotation procedure is to perform the integration 
over metrics of positive definite signature. However, unlike the case of fiat space 
quantum theory, this procedure is rather more than giving a mathematical mean- 
ing to the functional integral. The reason is that not all Lorentzian metrics have 
Fuclidean sections, and vice versa. Therefore, if the integral is performed over 
all Fuclideanized metrics then some Lorentzian metrics will have been omitted, 
whilst if the integration is taken over all Fuclidean metrics then this will include 
some which do not admit Lorentzian sections, and are thus not physical. Hawk- 
ing [3] has proposed a radical solution to this problem: the integral should be 
performed over all Fuclidean metrics, and this can be regarded as being rather 
like a contour integral in the space of all complex metrics in which the contour 
has been rotated from the Lorentzian to the Fuclidean section. Thus individ- 
ual metrics in the functional integral need not necessarily have a direct physical 
interpretation and only Z itself should be analytically continued back to the 
Lorentzian regime at the end. Consequently, the Fuclidean quantum-gravity 
path-integral is defined as 

Z = J D[g]e-'- , (21) 

M 

with the factor y/—g replaced by —iy^ and the Fuclidean action Ie = —il. 

But the problem is not solved yet because the Fuclidean gravitational action is 
not positive definite even for real positive definite metrics. Consider the situation 
of a compact manifold M without boundary, so that there is no surface term in 
the action (which is to be the case for the remainder of the article) and perform 
a conformal transformation gmn dmn = ^'^dmn under which 

i?^ i? = 0-2i?_60-3nO . (22) 
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Then the action becomes 



Ie[9] 



1 



/ 



mn 



(23) 



IQttG 



M 



One sees that /_£;[</] can be arbitrarily negative by choosing a rapidly varying 
conformal factor 0. It has been shown how this problem can be circumvented in 
the case of the functional integral over all metrics which approach flat Euclidean 
space at infinity (asymptotically Euclidean metrics). The integral (21) can be 
decomposed into an integral over all such metrics with Ricci scalar i? = 0, and 
an integral over the conformal deformations of these metrics. According to the 
positive action theorem the action of an asymptotically Euclidean metric with 
i? = is positive or zero, and vanishes if the metric is fiat [15]. The troublesome 
behavior has therefore been isolated in the conformal degrees of freedom, and 
by rotating to an integration over complex conformal factors, which makes the 
kinetic term for the conformal factor positive, this can be rendered convergent 
also [16]. 

A second serious problem is connected to the fact that the coupling constant 
has dimension of mass to the power (n — 2) and suggests that the theory is 
not perturbatively renormalizable above two dimensions. 



There are attractive physical reasons for considering 4-geometries with different 
topologies. That is to say, there can occur quantum fiuctuations of the metric not 
only within a given topology but from one topology to another suggesting a foam- 
like structure of spacetime on the scale of the Planck length [3]. This structure 
might give rise to an observable effect like the gravitational decay of baryons 
or muons. This can be caused, when they fall into gravitational instantons or 
virtual black holes and come out again as other species of particles. Evoked by 
this spacetime foam one could ask for quantum transition amplitudes between 
states specified by disconnected as well as connected 3-manifolds and multiply 
connected as well as simply connected ones. Unitarity would then suggest that 
the Euclidean functional integral contain a sum over topologically nontrivial 4- 
manifolds into which these 3-manifolds can be embedded. 

One can imagine explicitly implementing this sum by first making a list of 
all physically distinct 4-manifolds, which are not diffeomorphic to each other. 
But it can be shown, that there is no classification scheme for n-topologies, 
n > 4, because there is no universal algorithm to decide when two entries on 
the list are the same manifold [17]. Subclasses are classifiable: if one imposes 
the condition that the manifolds be simply connected and admit a spin structure 
(vanishing 2"''^ Stiefel- Whitney class) then they can presumably be classified (up 
to homotopy) by two topological invariants, the Euler characteristic x s-nd the 
Hirzebruch signature r. In the case of differentiable Riemannian manifolds, x and 
T can be expressed in terms of integrals involving the metric and the Riemannian 



2.2 Topology 



13 



curvature tensor as 



X = TT^T^ / ^'^v^^a6cdi?«™.e"'''e^'"'^ (24) 
T = ^ / d\^gK,,,W'^^e^'-^-^ . (25) 

In a recent paper a wider class of manifolds, called conifolds were shown to be 
algorithmically decidable in four dimensions [18]. 

2.3 Measure 

Several proposals for the nonperturbative integration measure have been sug- 
gested for quantum gravity. All of them are based on formal arguments which 
explain the variety of obtained results. Here, the origins of only a few are out- 
lined, without any pretence of completeness. 

The definition of the distance between two nearby geometries, whose differ- 
ence Sgij is invariant under reparametrizations and ultralocal, is [6] 

= j d\d)j6g,,{x)G'^"{x,y)6gu{y) , (26) 
where G^-'^'{x,y) is the local bitensor density 

G'^"{x, y) = v^(5V' + 5 V + 79''9"n^ - y) (27) 

and 7 7^ I but otherwise arbitrary. The determinant of the infinite continuous 
matrix G^^^'{x, y) can be computed at a formal level to get for the measure 

2 {n — 4:){n-\-l ) 

D[g] = n Gix)^ n dgtjix) = const ]^5'(a;) » ddtjix) , (28) 

X i<j X i<j 

which in four dimensions, n = 4, is simply a constant. 

Other results have been obtained demanding scale-invariance [19] and also in 
the Hamiltonian approach [20] 

D[g] = I[9(^)~'^I[d9^A^) ■ (29) 

X i<J 

To summarize, the discrepancies between various measures seem to lie in the 
definition of the product JJ^ leading to different powers of in the pre-factor. 
Although a rigorous definition of the measure on the space of geometries has 
been given allowing for an infinity of true Lebesgue measures the problem of 
which is the correct measure for quantum gravity is not yet solved at all [21]. 
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3 Simplicial Quantum Gravity 



Faced with these problems, various methods have been proposed to compute 
quantum-gravity path-integrals approximately. A well known way of doing so 
is to use simplicial manifolds leading to a lattice formulation of gravity called 
Regge calculus [10]. The Regge calculus needs a skeleton-like structure for the 
spacetime manifold and describes its geometry without using coordinates. 

3.1 Simplicial manifolds 

It is useful to begin by summarizing certain definitions and theorems on simplicial 
manifolds [18, 22]. 

Let Ui, ^2 . . . , Vn+i be affinely independent points^ in R""'"^. An n-simplex 
is the convex hull of these points: 

n+1 n+1 

= {x\x = J2 ^^v^; > 0; E = 1} • (30) 

8 = 1 8 = 1 

A 0-simplex is a point, a 1-simplex is a line segment or edge, a 2-simplex is a 
triangle, a 3-simplex is a tetrahedron. Higher dimensional simplices are general- 
izations of tetrahedra to higher dimensions. A simplex spanned by a subset of 
the vertices of an n-simplex is called a face. A simplex is uniquely determined by 
its vertices. This property is very important for both using and understanding 
simplicial complexes. 

A simplicial complex K is a topological space \K\ and a finite collection of 
simplices s'^ such that 

1. I I is a closed subset of some finite dimensional Euclidean space, 

2. if ^ K , then all faces of s'^ belong to K , 

3. if s'^ , G K , then either s'^ r\ = ^ ov s'^ r\ is a common face of both s'^ 
and . 

The topological space \K\ is the union of all simplices in K. The dimension of 
a simplicial complex is the largest dimension of any simplex contained in the 
complex. 

Thus, a simplicial complex describes both the building blocks of the space \K\ 
and gives the rules for how these building blocks are connected. Consequently the 
simplicial complex completely describes the topology of the space \K\. Finally, as 
each simplex in K is uniquely determined by its vertices, the simplicial complex 
itself is uniquely determined by the vertices and the connection rules. It is clear 
that this property is especially valuable for computational purposes. 

Simplicial complexes can describe topological spaces containing subspaces of 
different dimension, compact and noncompact spaces, and spaces with boundary. 

"'^A set of points {t;,} is afEnely independent if the vectors {vi — Vj^, Vi ^ vj, are linearly 
independent for all Vj . 
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Thus, there are simphcial analogs ol various standard definitions in topology. 
A pure simplicial complex is one in which every lower dimensional simplex is 
contained in at least one n-simplex where n is the dimension of the simplicial 
complex. A compact simplicial complex is one which contains a finite number 
of simplices. A connected simplicial complex is one in which any two vertices 
are connected by a sequence of edges. A nonbranching simplicial complex is a 
simplicial complex of dimension n in which every (n — l)-simplex is contained in 
at most two n-simplices. 

At this point it is possible to define special sets of pure nonbranching simpli- 
cial complexes. They are given the name of pseudomanifolds. A pseudomanifold 

is a pure nonbranching simplicial complex such that any two n-simplices can 
be connected by a sequence of n-simplices, each intersecting along some (n — 1)- 
simplex. 

In order to study subsets of pseudomanifolds, two more definitions are re- 
quired. These definitions are used to characterize the local topology of simpli- 
cial complexes and thus provide the means for defining simplicial equivalents of 
smooth manifolds. The star St[v) of a vertex v is the complex consisting of all 
simplices that contain v. The link L[v) of a vertex is the subset of the star of v 
consisting of all simplices in the star that do not intersect v itself, cf. Figure 2. 

Now the simplicial equivalents of smooth n-manifolds can be defined. A 
simplicial n-manifold is an n-pseudomanifold for which the link of every vertex 




Figure 2: The star of a given vertex v consists of the interior, edges, and vertices 
of those triangles which intersect v. The link of v is composed of the heavily 
drawn edges and their vertices in the Figure. It is topologically equivalent to . 
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is a simplicial (n — l)-sphere. A necessary and sufficient condition for the link 
to be a simplicial (n — l)-splLere is for the star of a vertex to be a simplicial n- 
ball. The main reason for defining simplicial manifolds in the above way is that 
it is a homogeneous definition, namely, no vertex of the simplicial complex has 
preferred treatment and the links of the vertices are all homeomorphic. This is 
similar to the idea of a topological manifold where each point has a neighborhood 
homeomorphic to a ball. 

In order to discuss the connection of smooth spaces to simplicial spaces, the 
first concept needed is that of a triangulation. A triangulation of a manifold Af" 
consists of a simplicial manifold and a homeomorphism t : \K^\ — > Af". Any 
topological space that admits a triangulation is homeomorphic to a polyhedron. 
Therefore, spaces allowing a triangulation are nice in the sense that they have 
the same properties as polyhedra. The most straightforward way to triangulate 
an n-dimensional manifold is to divide it into a sufficient number of hypercubes. 
Then draw the set of all diagonals (face, body, hyperbody diagonals) from each 
vertex to the other 2" — 1 points in the hypercube, i.e. the degrees of freedom per 
point are 2" — 1. Consequently the hypercube consists of n\ identical n-simplices 
and the coordination numbers are equal for each vertex. It is remarkable that in 
spite of the high symmetry the coordination numbers N^/Nn for the hypercubic 
triangulation are close to those of a random triangulation [23], which in general 
is desirable to use. Here Ni denotes the number of z-simplices, i = 0, . . . , n. 
A 2-dimensional example of this so-called regular triangulation is illustrated in 
Figure 3. 

In four dimensions we have 24 4-simplices per hypercube [24], with A^i = 
15iVo, N2 = 50iVo, N3 = 36iVo, and N4 = 24iVo. Thus, we have 15 edges k 
per vertex, which can be labeled in the following way. Taking the initial point 
of some 15 edges as the origin of a Cartesian coordinate system, we assign the 
coordinates (^4, ns, ^2,^1), rii G {0, 1}, of the endpoint to each edge and indicate 
li such that i = rii -\- 2n2 + 2^n3 + 2"^n4, cf. Table 1. 



12 3 1 




12 3 1 



Figure 3: A regular triangulated 2-torus is represented as a rectangle with 
opposite sides identified. 
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edges 


face diagonals 


body diagonals 


hyperbody diagonals 


1 (0,0,0,1) 


3 (0,0,1,1) 


7 (0,1,1,1) 


15 (1,1,1,1) 


2 (0,0,1,0) 


5 (0,1,0,1) 


11 (1,0,1,1) 




4 (0,1,0,0) 


6 (0,1,1,0) 


13 (1,1,0,1) 




8 (1,0,0,0) 


9 (1,0,0,1) 


14 (1,1,1,0) 






10 (1,0,1,0) 








12 (1,1,0,0) 







Table 1: Notation for labeling edges within a hypercube. 



3.2 Regge action 

To translate the path integral (21) into a concrete sum over simplicial histories 
it is necessary to be able to associate a metric and action with any K^. Thus 
in order to proceed, metric information must be attached. The easiest way 
to do this is to require that the metric on the interior of each n-simplex is 
the Euclidean metric, that is all n-simplices in the simplicial complex are flat. 
With this requirement, the geometry of each n-simplex is completely fixed by 
specifying the lengths li of all of its edges. Not every assignment of edge lengths 
is consistent with the simplices having fiat interiors. The triangle inequalities 
and their analogs for tetrahedra and 4-simplices must be satisfied. In more 
mathematical terminology, the simplicial complex becomes a piecewise linear 
manifold, that in general does no longer fit in R". It follows that the geometry 
of the simplicial complex is also completely fixed by specifying the lengths of 
all edges li in the complex and one anticipates that all geometrical quantities 
such as volume and curvature can be expressed completely in terms of the edge 
lengths. 

Indeed this is the case. Consider an n-simplex and define the n vectors 
ei which start with the vertex and proceed to the vertex z, i.e. the vectors 
ei, . . . , e„ span the n-simplex. Its volume Vn can be expressed by the generaliza- 
tion of Tartaglia's formula for a tetrahedron 

K = ^\/det{e, ■ ej) , (31) 

where the n X n-matrix of scalar products may be written in terms of the quadratic 
edge lengths qij = \ei — e^p between vertices i and j by 

• fij = ^{qo, + qoj - q,j) . (32) 

Somewhat less obviously, curvature can also be expressed in terms of the edge 
lengths. As the metric on the interior of each n-simplex is fiat, it is clear that 
the curvature of the combinatorial space is not carried on the interiors of the 
n-simplices. Rather, it turns out to be concentrated on the (n — 2)-simplices of 
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the simplicial complex. This becomes directly apparent in two dimensions, in 
which curvature is concentrated on vertices. The curvature is defined locally by 
the amount of rotation of a vector parallel-transported around an infinitesimally 
small loop. Since the space is fiat in each triangle, the rotation of the vector 
is zero unless the loop contains a vertex v inside. Then the circum-transported 
vector rotates by an angle S[v). Such a region where the (5-function-like curvature 
resides is called a hinge, and the rotation angle S[v) is called deficit angle. Let 
m denote the number of triangles in St[v)^ then the deficit angle associated with 
this vertex is given by (cf. Figure 4) 

m 

6{v) = 27r-J2e, , (33) 

8 = 1 

and 9i is the angle between two unit vectors that lie in adjacent edges of the i^^ 
triangle. Indeed the sum of the deficit angles over all vertices in the 2-manifold 
Af^ yields the Euler characteristic 

Y: S{v) = 27rx{M') , (34) 

as required by the Gauss-Bonnet theorem. 

The above discussion can be easily generalized for simplicial manifolds with 
dimension n. In this case the hinges are (n — 2)-simplices. If the index i = 
l,...,m sequentially labels adjacent n-simplices in S't(s"~^), then the deficit 
angle is given by 

m 

6{s--') = 27r -J2e, , (35) 

8 = 1 




Figure 4: A 2-dimensional simplicial geometry is a net of fiat triangles together 
with an assignment of lengths to their edges. The curvature is concentrated at 
the vertices and is measured by the deficit angle 6[v) which a parallel-transported 
vector experiences. 
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where 6i is now the dihedral angle constructed from two unit vectors normal to 
that lie in the adjacent (n — l)-simplices of the i^^ n-simplex. The dihedral 
angle can be computed in terms of the edge lengths by elementary trigonometry 
in any dimension. 

Finally, given the above definitions, the Regge action for Einstein gravity 
with cosmological constant A for a closed pure nonbranching complex is 
[10, 11, 25] 

lR = -m'p J2 26{s^-')V{s^-') + ml J2 '^^^i^'") ^ (36) 

where the first sum is over all (n — 2)-simplices in the simplicial complex and 
the second is over all n-simplices in the complex. is the n-volume of the 

indicated n-simplex and mp = [l/lGirG]^ is the bare Planck mass. The Regge 
action for compact simplicial manifolds with boundary can also be formulated 
entirely in terms of edge lengths. Essentially one adds the appropriate discretized 
form of the boundary term (18) that appears in the continuum action [11]. A 
detailed analysis shows that in the classical limit the Regge action converges 
towards the Einstein-Hilbert action if the fatness (defined in the next Chapter) 
of each simplex is finite [26]. The classical continuum limit is reached with 
increasing the number of vertices A^o — ^ and simultaneously decreasing the 
edge lengths qi — )• 0, i.e. the local lattice spacing becomes smaller than the local 
radius of curvature. 

3.3 Group action for simplicial gravity 

Given a simplicial lattice, to each edge of its dual lattice a Poincare transfor- 
mation can be assigned. With the new variables from these transformations, it 
is possible to build an action that reduces to the Regge action Ip in the small 
curvature limit [27]. 

There is a general procedure for the construction of a dual lattice, whose cells 
are Voronoi polyhedra [23, 25]. The Voronoi polyhedron dual to vertex P is the 
set of all points that are closer to P than to any other vertex. In this way a 
A;-simplex in the lattice is dual to an (n — A;)-polyhedron in the dual lattice. 

Let us consider a particular hinge with its vertices Pi,...,P„_i and 

its star St[s^~'^) consisting of the set of n-simplices s^, . . . , sj^ characterized by 

= {Pi, . . . , P„_i, Qi-i^i^ Qi.i+i}, {i = 1, • • • , rn)- Eurther let us choose in each 
n-simplex a Lorentzian frame and denote the coordinates of the vertices of 
in this frame by 

P, = {j/»(z), a = l,...,n} 
Q^-l,^ = a = 1, . . . , n} , (37) 

and also for each vertex of the dual lattice Di = {a;"(z), a = l,...,n} (cf. 
Eigure 5). Eurthermore, we associate to each edge of the dual lattice a Poincare 
transformation 

U{i, i + l) = {U^{i, I + 1), [/"(z, z + 1)} , (38) 
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1 



Figure 5: Part of a 2- dimensional lattice and its dual. 



and define [/(z, z + 1) by demanding 

U^{t, I + l)y]{i + f ) + z + f ) = j/;(z) (39) 

and analogous for the other vertex coordinates. The arbitrariness of the choice 
of the reference frame in each n-simplex becomes in the dual lattice the gauge 
invariance under local Poincare transformations n(z): 

[/(z,z + i) ^ n(z)[/(z,z + i)n(z + 1) 

x^{i) Jll{i)x\i) + Jr{i) . (40) 

By suitable fixing the translational part of the gauge it is always possible to set 
x°'{i) = Vz. Therefore, the plaquette variable on the dual lattice around the 
hinge 

W,{s^-^) = U{i, z + l)U{i + 1, z + 2) . . . [/(z - 1, z) (41) 

leaves the coordinates of all vertices P,- of unchanged and becomes a rotation 
by the deficit angle 6{s^~'^) in a plane perpendicular to 

This leads to the following proposal for a "compact" gravitational action on 
the lattice in terms of plaquette variables 

n-2 

gn-2 j = l 

= -m^ ^ sin,5(s"-2)y(s"-2) . (42) 

gn-2 
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It is obvious that (42) reduces to (36) in the hmit of small deficit angles, further 
it bears a striking resemblance with the continuum action written in terms of 
differential forms. 

In deriving (42) the Levi-Civita connection has been implicitely used for the 
parallel transport within a loop. But Ic can also be defined for general connec- 
tions, including ones with nonvanishing torsion and even ones not determined by 
the metric [28]. 

One disadvantage of the compact action is, that it does not have the property 
Ir{M"') = Ir{M^) + /i^(Af^) like the Regge action, whereby we have defined 
Af" = + M2 as the manifold obtained by joining together the manifolds 
and M2 along M^~^ . It is therefore not clear, whether it gives rise to a unitary 
theory in the form of physical- or refiection positivity [28]. 

3.4 Simplicial path integral 

Given a method of associating a geometry and an action with any complex 
the Regge equivalent of the path integral (21) is 



where the integration over all n-geometries without boundary has been re- 
placed by a sum over distinct simplicial complexes and integrating over their 
quadratic edge lengths q. Usually one abandons summing over different but 
contents oneself with only one fixed toplogy. Anyhow, arguments have been 
pointed out that a sum over topologies cannot give a finite functional integral, 
because the number of complexes within a given topology increases too rapidly 
as a function of the number of simplices [28]. 

Two somehow complementary methods have been established to evaluate the 
simplicial path integral. 

a) The Regge approach: The triangulation is kept fixed and the edge lengths, 

which play the role of infinitesimal geodesies in the continuum are allowed 
to vary. The path integral (43) reduces to a summation over different 
configurations [q] [II, 25, 29, 30]. 

b) The dynamical triangulation approach: This method represents an al- 

ternative and complementary approach to what will be discussed here [31]. 
The simplicial manifold consisting of a set of equilateral simplices is locally 
changed by applying elementary moves (e.g. the well-known Alexander 
moves [32]) on the incidence matrices, i.e. one performs a summation over 
all possible triangulations of the manifold. In two dimensions this method 
is equivalent to matrix models, but in four dimensions problems seem to 
appear, because it is not possible to triangulate fiat space with equilateral 
4-simplices and therefore to reach the classical continuum limit directly. 




(43) 
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Since within the Regge approach the squares of the edge lengths, are hnearly 
related to the corresponding measure to D[g] would be D[q] = Hg; dfJ-i^i) 

indicating the freedom available in choosing the measure on the space of edge 
lengths. Investigations concerning regularly triangulated hypercubic lattices sug- 
gest [24, 25, 26, 29, 30] 

D[q]=I[<ir'dqiHq} , (44) 

with = 1 if the generalized triangle inequalities are fulfilled and = oth- 
erwise. The parameter cr > determines the behavior of the measure under 
rescaling. The value (7=1 corresponds to the uniform measure (28), whereas 
for cr = a scale-invariant measure results (29). 

The functional integral over metrics has now been reduced to the product of 
integrals over edge lengths as the metric is discretized. This has the consequence 
that diffeomorphism invariance is destroyed. The recovery of the spacetime dif- 
feomorphism group is of vital significance in discussing the continuum limit and 
can be approached from several perspectives. For example, a triangulation of a 
manifold M can be viewed continuous map r from an abstract simplicial 

complex K to M. Then each diffeomorphism ^ of Af generates a new triangula- 
tion (po T : K ^ M . Unfortunately, the discretized Regge action is not invariant 
under the change in edge lengths, there remains only an approximate invariance 
group, which becomes increasingly accurate with increasing refinement of the 
lattice. This is revealed by the fact that the Hessian matrix constructed out of 
second derivatives of the action with respect to edge lengths is generally nonsin- 
gular [11] (unlike the continuum counterpart). Whether or not such a structure 
would suffice for quantization is not known. But in any event, there is a clear 
danger of overcounting the metric modes and thus a type of gauge fixing for 
the triangulation was developed [33]. Given a complex to each metric on M 
corresponds just one triangulation r obtained by minimizing the deviations of 
the edge lengths from their mean. This is realized by dynamical triangulation, 
which is plagued by other problems like ergodicity of the path integral [34]. 
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4 Entropy- dominated Phase 



Different formuiations of tlie simpiiciai action iead to tlie same ciassicai con- 
tinuum iimit, but it is not ciear to wliat extent tliey differ after quantization. 
Unfortunateiy, tlie computation of tlie path integral is not straightforward, be- 
cause the action is still unbounded and the Regge calculus admits no simple 
means of dealing with conformal modes and hence controlling this deficiency. 
Indeed, both the original Regge action (36) in four dimensions, n = 4, 

/« = -/3^AA + A^V; (45) 

t s 

as the action (42) with the compact deficit angle contribution and an additional 
cosmological term 

Ic = -/3j2Atsin6t + Xj2Vs (46) 

t s 

can be made arbitrarily negative if some of the 4-simplices have near-zero 4- 
volume Vs but contain very-large-area triangles Aj, whose deficit angles are pos- 
itive. It is reasonable to set an appropriate cut-off via a lower limit / on the 
fatness 

<Ps = C' C = 24, (47) 

of each 4-simplex, ^ / > [26]. The fatness is maximal for equilateral 
simplices and goes to zero for collapsing ones. 

But even without any cut-off the unboundedness of the action need not to 
render the path integral ill-defined. The reason is that entropy effects coming 
from the measure could cure this pathology [29]. This can be seen if the path 
integral is rewritten as 

oo 

Z = j dl n{I)e-^ , (48) 

— oo 

where n[I) gives the density of states for the action /. If n[I) goes faster to zero 
with / — )• — oo than the exponential diverges, then Z would exist. Numerical 
computations indicate that the entropy of the system can indeed compensate the 
unbounded action [25, 29, 30]. 



4.1 Monte Carlo method 

A possible nonperturbative method to approximate expectation values of observ- 
ables O 

(O) = ^ f UMqi)Om)e-'^^'^^^ (49) 

is by using numerical Monte Carlo (MC) techniques that do not rely on an 
expansion in a small parameter. A finite discrete system is considered, which 
is no serious restriction from the viewpoint of computer simulations, because a 
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continuous variable is anyhow discretized by the computer's digital accuracy. For 
a finite system, (49) is replaced by 

{0) = ^Y.Oke-'' ^ (50) 

k 

where the possible configurations {qi}k are labeled by integers k and therefore 
Ik = I{{(li]k)- The division by the partition function Z is to normalize expec- 
tation values, such that (1) = 1. Since the total number of configurations is 
exorbitantly large one performs an importance sampling, so that configurations 
are generated with probability 

P{{qi]k) = const 6-^'= . (51) 

The constant is determined by the normalization condition Y,kP{{^i}k) = 1- 
Configurations with (51) as their equilibrium distribution can be produced by 
a Markov process. The elements of the Markov chain are the configurations. 
They are sequentially generated such that {qi}j depends only on the previous 
configuration {qi}j-i. The transition probability to create {qi}j from {qi}i is 
given by the matrix element Wij = W{i — > j) > 0, which is required to satisfy 
the following properties: 

i) Ergodicity: 

e-^' > , e"^^ > (52) 
implies the existence of an integer number m > 0, such that {W"^)ij > 0. 

ii) Normalization: 

E = 1 . (53) 

i 

iii) Balance: 

J2w,,e-'' = e-'^ (54) 

i 

means, the Boltzmann weight factors form eigenvectors of the matrix W. 

An ensemble is a collection of configurations such that to each configuration a 
probability Pi is assigned and Pi = 1. The Boltzmann ensemble is defined 
to be the ensemble with probability distribution 

Pf = const e-^' , Y^Pf = 1 . (55) 

i 

An equilibrium ensemble Ee of the Markov process is defined by satisfying 
WPe = Pe- The central point for the MC method is, that under the condi- 
tions i), ii), and iii) the Boltzmann ensemble is the only equilibrium ensemble of 
the Markov process (see [29] for a proof). There are many ways to construct a 
Markov process satisfying these conditions. In practise most MC algorithms are 
based on the condition of detailed balance: 

W,,e-^' = Wj,e-^^ . (56) 
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Using the normalization condition ii) = 1? detailed balance immediately 

implies balance iii). A popular choice is the Metropolis algorithm [35] because 
of its computational simplicity. This procedure is described below for the Regge 
skeleton. 

A new configuration with Ik is proposed by varying the edge lengths of the 
lattice by a small amount according to the measure (44). If it is not in the 
Euclidean sector, i.e. the triangle inequalities are not fulfilled, or if the fatness 
(ps is lower than the limit /, it is rejected. Otherwise one has to check the 
Metropolis condition: the new configuration Ik is always accepted for Ik < Ij 
and for Ik > Ij it is accepted with probability exp(/j — Ik). 

A MC iteration is finished, when a new value is proposed for every edge in 
the lattice. In the limit of an infinite number of iterations the expectation values 
are given by 

{0)=li^^j:Oki{q^}). (57) 

^ k=l 

In practice, the generated configurations are averaged over a finite number of 
iterations after thermalization to approximate the considered expectation value. 

4.2 Results 

MC simulations have been performed to investigate the entropy-dominated phase 
for both proposed simplicial gravitational actions, (45) and (46), with the fol- 
lowing features and parameters. 

Lattice: The regularly hypercubic triangulated simplicial manifold is fixed to be 
homeomorphic to a 4-torus T"^. The total number of vertices is A^o = 3"^ X 8. 

Fatness: To facilitate the simulations a lower limit / = 10"'* is set for the fatness 
(ps of each 4-simplex. This scale invariant cut-off reduces the configuration 
space and thus allows to reach an equilibrium after a reasonable number of 
MC iterations. Previous computations suggest that such a restriction does 
not affect the behavior of the system in the well-defined phase [30]. 

Measure: Investigations of a 1-parameter family of measures (44) indicate the 
stability of the entropy-dominated region against variations of the measure 
[30]. As the simplest choice the uniform measure with cr = 1 is used, but 
universality of the results is expected. 

Couplings: The bare cosmological constant appearing in the gravitational ac- 
tion (45) or (46) is set to A = 1 by fixing the average lattice volume. /3 
determines the gravitational coupling and is varied in the range [0.07, 0.15]. 

Iterations: 100k - 300k MC sweeps inclusive thermalization have been applied 
for each value of /3. 
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Examining the scaling properties of the action and the measure, a very useful 
identity for checking the accuracy of numerical computations is obtained [25]. 
The path integral 



obeys 



which implies 



Z(/3,A) 

2A(^K)-/3(^AA) = ^TiVi 



(58) 

(59) 
(60) 



an equation easily to verify in simulations. The correspondence between analyt- 
ical and numerical results is shown in Figure 6 for the Regge action. The same 
behavior is found for the compact action Ic- 

A quantity of physical interest is the average curvature 



n 



IV9 



(61) 



which can be understood as effective cosmological constant [30]. The variables 
At, y<;, q are expressed in units of the bare Planck length. 

The behavior of the average curvature as a function of the coupling /3 is de- 
picted in Figure 7. As /3 is varied a phase transition is found separating a smooth 



48A(K) - 50(3{ASt) 



15 



1 1 

-■o 


1 1 1 

o- ■«>o 


1 1 

o 


1 1 

o- 


1 1 


1 1 1 


1 1 


1 1 



0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 

/3 



Figure 6: For the 3"^ X 8 lattice with A = cr = 1 one has = 15, which is 
in perfect agreement with our simulation data, here exemplified for the Regge 
action Ir. Frror bars are smaller than the symbol size. 
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Figure 7: Average curvature using the Regge and the compact action Ir and /c, 
respectively, on a 3"^ X 8 lattice with A = cr = 1 for different gravitational couplings 
/3. In the well-defined phase most data points from both actions overlap. Error 
bars are smaller than the symbol size. 



from a rough phase of gravity. For sufficiently small /3 the curvature is small and 
negative and enters a region of large curvature for /3 > jScrit- How far the location 
of the transition point depends on the choice of the fatness, the measure, and 
the underlying triangulation is discussed in [30]. While the transition is rather 
discontinuous for the Regge action Ir^ it is smooth for /c, which might be im- 
portant for the continuum limit. This suggests that the appearance of sinSf in 
the action provides a dynamical mechanism to force the curvature to be small. 

A similar effect was reported when higher-order terms in the curvature, e.g. 

terms, are added to the Regge action Ir and simulated within the uniform 
measure. It was stated that the phase transition to the ill-defined region suggests 
some sort of multicritical behavior with a line of I**-order transitions without any 
higher-order term and a line of 2"''-order transitions with terms [25]. Unfor- 
tunately, this advantage for a continuum limit is accompanied by a violation 
of unitarity. Originally, the Finstein Lagrangian was replaced with one that 
involved squares of the Riemannian curvature because of perturbative renormal- 
izability. But this is achieved at the expense of acquiring a massive spin-2 ghost 
partner (Poltergeist) for the physical massless graviton [36]. 

The transition to the region of large curvature is accompanied by the tendency 
of the simplices to collapse into degenerate configurations. This is seen if one 
examines the expectation value of the fatness {(f>s)- Figure 8 shows {(ps) as a 
function of the coupling /3 for both actions. It turns out that in the well-defined 
phase the average fatness is much larger than the lower limit / = 10"'*. With 
the transition to the ill-defined phase the fatness decreases (suddenly) signaling 
the presence of collapsed 4-simplices (spikes) and a crumpled lattice with fractal 
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Figure 8: Fatness (^s) versus /3 in the 3"^ X 8 lattice with A = cr = 1 for the Regge 
action and the compact action. The cut-ofF for the fatness was set to / = 10"'* 
to allow equilibrium in the ill-defined phase. 



dimension less than four. This region of the phase diagram can only be observed 
for finite /, because without any cut-off the action diverges and simulations 
would not reach an equilibrium. 

The development of spikes is refiected in the formation of an inhomogeneous 
distribution of the edge lengths. Figure 9 depicts the histograms for a configu- 
ration of edge lengths below and above the transition point using the compact 
action Ic (46). While a homogeneous distribution is typical for the well-defined 
phase, it is rather inhomogeneous in the ill-defined regime. A similar behavior 
for the Regge action Ir (45) was found in [30]. 

To summarize, it seems that the system prefers configurations with distribu- 
tions of edge lengths in the well-defined phase leading to a small negative average 
curvature. Increasing the coupling /3 one reaches a critical value entering into 
the ill-defined phase characterized by degenerate configurations with collapsing 
simplices and very large positive curvatures because of the unbounded action. 
Note that this collapsed phase is the coupling region of the weak-field expansion 

In principle it could be possible that the well-defined region of the path 
integral is only a numerical artefact due to metastable states. Numerical sim- 
ulations examining inhomogeneous start configurations with very large positive 
curvature, however, indicate independence from the initial conditions, i.e. the 
average curvature returns to small and negative values for coupling parameters 
of the well-defined phase [30]. 

Another important task is to study the infiuence of the underlying triangu- 
lation on the entropy-dominated phase. The most natural setting would be a 
random lattice reproducing spacetime isotropy stochastically. Since random lat- 
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Figure 9: Histograms of the squared edge lengths q in units of the average squared 
edge length q = -^Y^iqi oi a, typical configuration below (upper picture) and 
above (lower picture) the phase transition applying the compact action Ic- 

tice simulations are technically cumbersome, most numerical computations have 
been performed on the regular hypercubic triangulation of the 4-torus. Never- 
theless, to investigate the effect of the lattice structure one can apply bary centric 
subdivisions and study such irregularly triangulated systems as the number of 
randomly inserted vertices increases. Recent computations reveal that the form 
of the incidence matrix and the local coordination number have an important 
infiuence on the critical coupling [37]. 
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5 Two-point Functions in Simplicial Quantum 
Gravity 



The previous discussion has dealt exclusively with averages of local operators, 
which provide a considerable amount of information. To understand the interac- 
tion mechanism of a theory one has to study correlation functions, which yield 
the mass spectrum of the system. This can be used to decide whether a massless 
graviton does exist or not. 



5.1 Weak-field limit 

For a 4-dimensional manifold, regularly triangulated as described in Section 3.1, 
the variation of the action (45) (A = 0) with respect to the edge length / is 

(9/4 

SIr = SiJ2 AA) = Y,SAS, = Y, -T^Si St , (62) 

t t t 

using the fact that At SSf = [10]. Hence the equations of motion are 

being the lattice analog of the Einstein equations. To derive the free two-point 
function the second variation of Ir around flat space, where the deficit angles all 
vanish, yields 

t 
t 

= i:(i:f«)(i:f«,. ,64) 

Now let (5i (z = 1, . . . , 15) be the amount of fiuctuation of U 

6k = , (65) 

then one can write (64) in the bilinear form 

S'^Ij^ = S^MS , (66) 

where (5 is a column vector with 15 components per lattice point and the matrix 
M corresponds to the inverse propagator. 

Expanding the fiuctuation at a point with coordinates (a, b, c, d) into 

periodic modes uji = exp(^) (i = 1,2,4,8) 

^(a,b,c4) ^ ^^^y^^^ f^^^y^^^ys, , (67) 
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one can write (66) in momentum space 

= 8lM:j8, . (68) 

By a unimodular similarity transformation 8^8' one can block diagonalize Af^^. 
It turns out that five modes are completely decoupled from the other ten and do 
not affect the dynamics of the system. Therefore, one considers only those parts 
of the vector 8\ whose degrees of freedom have the same multiplicity of ten as 
gran in the continuum theory. 

By means of a further transformation 8' ^ h the dynamical part of (68) is 

]-LR = h^2 J2 (1-cos— )f2^ ^^h-h^C^Ch. (69) 



i? is a 4 X 4 matrix, Iq the 6-dimensional identity matrix and C = C(cu) a 4 X 10 
matrix. For a concrete expression see [24]. 

In the continuum theory the Langrangian of the weak gravitational field hki 
with gauge breaking term Cm can be written as [40] 

C = -\cUuVumnrUmn + , (70) 

where 

V.lmn = \8k-ra8ln-\8H8mn (71) 

Cm = dnhmn " -Omhun • (72) 

To compare this expression with the result in the simplicial case, trace reversed 
variables are first employed 

hmn = hmn " -8mnhu , (73) 

and then the fluctuation variables relabeled such as 

hm+n = hmn , TTl ^ U . (74) 

Translating the variables into momentum space one finds 

L = hmh^ ^ kmh - h^C^Ch , (75) 

where km = idm is the momentum and G = G{k). By setting ki = — in the 
long wave length limit one can expand 

io, = l- ik + , (76) 

and gets 

2 (I -cos-) ^ ^ ^^ C ^ C. (77) 

8 = 1,2,4,8 ^« 8 = 1,2,4,8 

Inserting this in Equation (75) and comparing with (69) the continuum and 
lattice correspondence is clear (up to a trivial normalization factor) in the case 
of a weak gravitational field. 
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5.2 Correlation functions 



The continuum weak-field propagator can be formulated in terms of spin-2 and 
spin-0 projection operators [41] 

r p(2) _9n(0) 



with 



,(2) _ 1 (^PkPlf^ , PmP 
1 



^ ^{PkPm^ln + PkPnSlm + PlPmhn + PlPnhm) + 

1 12 

+ 77(4m'^te + hn^lm) " -^kl^mn + TT^PhPlPmPn , (79) 



p(o) _ _1 f -u 



^"f " 4f ] + -^klSmn + Tr^PkPlPmPn ■ (80) 



(v^i?(x)v^i?(j/)) ~ / / J4^i^(p,5)-^i^(p,5)^ , (82) 



klmn o \ 9 9 

O \ p p 

To probe the correlations in the spin-0 and spin-2 part separately one uses the 
fact that fluctuations in the volume are sensitive to the scalar channel 

(v^(x)v^(j/)) ~ - I d'p^e-^^(^-y^ , (81) 

whereas fluctuations in the curvature are due to presence of spin-2 particles 

p(2) . p(2) 

—-R{p,q)- 

[p — q) 

where R represents some momentum dependent vertex [42]. 

If A 7^ 0, the cosmological term implies a contribution to the propagator and 
the appearance of a tadpole diagram. With the gauge breaking term (72) the 
Langrangian becomes 

jC = \{1 + ^hu) + ^hklVklmnid^ - ^)hmn , (83) 

leading to almost the same propagator as before (78), except that p^ is replaced 
hy p^ -\- \ in the denominator. This propagator therefore corresponds to the 
exchange of a particle with mass m = \/\. At first instance the gravitational 
force becomes a Yukawa force with a range given by this mass. But it has 
been conjectured that by expanding around the correct solution of the classical 
equations of motion this mass term disappears, because in the presence of a 
cosmological term fiat space is no longer a solution of Einstein's equations. Up 
to now this has not been worked out even for simple cases [40] . 
In case of the Regge approach connected correlation functions 

G{d) = {0{x)0{y)), (84) 
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of a local operator O are measured at two points separated by the geodesic 
distance d = \x — y\ [42]. The geodesic distance d between two points x and y in 
a discrete space is the length of the shortest path joining them. 

Since the edge lengths are the dynamical variables of the theory a logical 
analog of (78) would be the 15 X 15 matrix 

Gi\d) = {qa{d)q,{0)), . (85) 

The indices a, b label the different types of edges within a hypercube. The squared 
edge lengths are used mainly for technical reasons. Some care must be taken to 
interpret (85) as a propagator. From perturbative calculations considering small 
fluctuations around flat space (cf. previous Section) it is already known that 
the free inverse propagator on the lattice has zero eigenvalues. One can argue 
that these zero modes are not present on a general curved lattice since diffeomor- 
phism invariance is lost. Moreover, the formulation of gauge transformations and 
gauge fixing is merely worked out for lattice gravity because the Regge calculus 
admits no natural means of introducing a ghost determinant or a gauge breaking 
term. But the fact that the lattice theory is already formulated in a coordinate 
invariant manner at least legitimates correlations of the type (85). Also correla- 
tions between deficit angles G\^{d) = {6a{d)6b[0))c^ with the indices a, b labeling 
the 50 triangles attached to a vertex within one hypercube, are confronted with 
the above problems, but can give insight into the interaction mechanism of the 
system. 

The simplicial analogs of the invariant two-point functions (81) and (82) 
correspond to the volume correlations 

G(^Hd) = {j:Vs E Vs')c (86) 

and the curvature correlations 

G'(«)(J) = (EMt E'^t'^t')c- (87) 

Here "^Zsdvo denotes the sum over all volumes of 4-simplices meeting at the 
vertex Vq and "^ZtDvo ^t^t collects the contributions of all triangles at Vq. The 
vertices Vq and Vd are separated by the geodesic distance d. 

5.3 Results 

MC simulations have been performed to compute several correlation functions 
for both simplicial gravitational actions, (45) and (46). In the presentations the 
index labeling the vertices in the long lattice direction is taken as a number for 
the actual distance, henceforth called index distance, between local operators. 
This seems to be reasonable because in the well-defined phase with its small 
average curvature the index distance is presumably a good approximation of the 
true geodesic distance. Further, it avoids a time consuming determination of 
geodesic distances via random walk on the fiuctuating simplicial lattice. 
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Lattice: A hypercubic triangulated 4-torus with A^o = 3"^ X 8 and A^o = 4"^ X 16 
vertices is used. Because of the inherent periodicity of the 4-torus the 
maximum index distance for two-point functions is 4 and 8, respectively. 

Fatness: The fatness cut-off has been set to / = fO""* and allows the study of 
correlations even for gravitational couplings /3 > /3c- 

Gravitational coupling: [3 is varied within the range [0.07,0.15] for the small 
lattice and [0.1,0.14] for the large. 

Measure and cosmological constant: They are set equal, cr = A, so that (Vs) 
takes the same value N1/2N4 for cr > at /3 = 0, cf. (60). In particular, 
two values of cr = A were chosen for every lattice size: 0.25 and I.O for the 
smaller and 0.1 and 1.0 for the larger. 

Iterations: lOOk - 300k and 50k MC sweeps were produced for the 3"^ x 8 lattice 
and the 4"^ X 16 lattice, respectively. The 3"^ X 8 lattice was first used to 
observe the tendencies of correlations between squared edge lengths (85) 
and between 4- volumes (86). 

Figure 10 exhibits the correlations ^ of squared edge lengths with edges parallel 
to the long lattice direction. For the Regge action in the well-defined phase, 
/3 < /3c ~ 0.1055, correlations are very short ranged and more or less zero for 
d > 2. Arround the critical point /3c edges are correlated with a long range and 
shallow slope. For /3 > /3c, correlations change considerably and tend to oscillate. 

The situation for the action Ic (lower plot of Figure 10) is not very different, 
but correlations have larger range for the corresponding /3 values and there is 
no zigzag behavior seen. The oscillating domain for Ic is possibly reached for 
/3 > 0.15, where the average fatness becomes small enough. 

In Figure II our results for the largest eigenvalue 7^^^ of the edge correlation 
matrix G^J (85) are drawn. Here all 15 edges are involved in contrast to the last 
plot dealing with edge lengths in the long direction. The maximum eigenvalues 
indicate two different types of correlations corresponding to the two phases of 
the system when using Ir^ whereas for the case of Ic no unequivocal separation 
into two kinds is possible. Oscillations become more pronounced with increasing 
gravitational coupling. One should remark that considering the correct geodesic 
distance could change all pictures. 

In constructing volume correlations, the sum over 4-volume elements at site 
V in (86) was restricted to one hypercube in the direction of d. This keeps the 
observable more localized and leaves the results practically unchanged. Volume 
correlations G^^^ are presented in Figure 12 and behave similar as the eigenval- 
ues 7^^^ of the edge length correlations. For the Regge action in the smooth 
phase they vanish for d > 2 and oscillate with increasing /3 in the rough phase. 
The volume correlations for the compact action are nearly proportional to the 
edge correlation eigenvalues, cf. Figure II. It is remarkable that the volume 
correlations for Ic are strictly positive in comparison to those for Ir. 
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Next, in Figures 13 - 15 for the Regge action in the well-defined phase the 
edge and volume correlations are contrasted with those obtained for a different 
choice of the cosmological constant, A = cr = 0.25. Since it was demonstrated in 
earlier work [30] that one-point expectation values do not change very much in the 
range 5 cr ^ 1, in this way one can directly study the effect of the cosmological 
constant A on the two-point functions. This subject will be discussed now on the 
larger lattice. 

In general, for correlations (84) one expects the functional form [42] 

G{d) = Cj^e--' , (88) 

which was investigated quantitatively on a 4"^ X 16 lattice. Figures 16 - 18 display 
the results of edge, volume and curvature correlations for the Regge action. In 
our analysis several fit philosophies have been tried. One can attempt to fit the 
correlations to a pure exponential function or try a power law. It turned out, 
that the exponential decay is more prefered suggesting finite effective masses. 

For the maximum eigenvalues of G^^^ the fit parameters are given in Table 2. 
In both cases, A = cr = 1.0 and A = cr = 0.1, a decrease of the effective mass 
m^') towards the critical coupling is visible, consistent with a massless particle at 
criticality. On the other hand, a decrease due to smaller cosmological constant 
is not observed. The same behavior turns out for the maximum eigenvalues of 
the deficit angle correlation matrix G^^^ . 

The volume correlations G^^^ show a decreasing mass m^^^ as A is reduced, 
cf. Table 3. When approaching the critical point for constant A = cr = 1.0 
or A = cr = 0.1, respectively, one finds contrary effects: m^^^ increases in the 
first case for /3 — )• /3c ~ 0.108, but decreases for /3 — )• /3c ~ 0.14 in the latter. 
Thus, infering a massless scalar particle as expected from weak-field calculations 
is hardly possible. 

The curvature correlations G^^^ in Figure 18 behave rather different com- 
pared to G^^^ or 7^^^. hi d = 1 all values of G^^^ are considerably below zero, 
which means that the local curvatures at neighboring vertices are strongly anti- 
correlated. Anyhow, for the fit procedure only index distances d > 2 are taken 
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Table 2: Fit parameters for the maximum eigenvalue of the edge length and 
deficit angle correlation matrices to a pure exponential. 
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Table 3: Fit parameters for the volume and curvature correlations to a pure 
exponential (C ~ 1). 

into account, because correlations at large distances should signal the occurence 
of massless gravitons. In the case A = cr = 0.1, m^^^ decreases similar as m^^^ 
and m^') when approaching the critical point. For A = cr = 1.0, G^^^ has strong 
fluctuations and it was impossible to fit correlations for /3 = 0.108. To get an 
idea of the local curvatures to be correlated in this case, see Figure 19, which at 
least gives some insight into the anticorrelated nature of the observable. It seems 
to be too daring to deduce the existence of massless spin-2 exchange particles 
from our data. 

To conclude, we have obtained first qualitative results for two-point func- 
tions of 4-dimensional simplicial quantum gravity [43]. Geodesic distances are 
not taken into account, which would be a considerable effort without changing 
things drastically in the well-defined phase. Although no convincing evidence for 
the particle spectrum predicted from weak-field approximations could be found, 
at least simulations with small cosmological constant and almost scale-invariant 
measure (A = a = 0.1) show a decrease of effective masses approaching the 
critical coupling. For computations applying the uniform measure and the cor- 
responding larger cosmological constant (A = cr = 1.0) statistics is not sufficient 
to yield clear signals for curvature correlations. The question of the infiuence of 
the cosmological constant on the graviton mass thus remains unanswered. 

Recently, interesting results were also reported for volume and curvature cor- 
relations with an additional term in the action [42]. There, an approximation 
for the geodesic distance was applied by measuring the relaxation length of a 
scalar field. These results are in favor of a massless graviton and thus qual- 
itatively different from our data presented in this Section. However, a direct 
comparison is not possible at this line because our notions of distances differ. 
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Figure 10: Correlations between edge lengths in the long lattice direction, , 
at different vertices separated by the index distance d for the Regge action Ir 
(upper picture) and for the compact action Ic (lower picture) with cr = A = 1.0 
on the 3"^ X 8 lattice. 
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Figure 11: Maximum eigenvalue 7^^^ of the edge length correlation matrix G^^jj 
for Ir (upper picture) and Ic (lower picture) with cr = A = 1.0 on the 3"^ X 8 
lattice. 
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Figure 12: Correlations between 4-volume elements for Ir (upper plot) and Ic 
(lower plot) with cr = A = 1.0 on the 3"^ X 8 lattice. 
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Figure 13: Correlation functions between edge lengths in the long lattice direction 
at different vertices separated by the index distance d. The gravitational coupling 
lies in the well-defined phase using the Regge action with A = cr = 1.0 (upper 
plot) and A = cr = 0.25 (lower plot). 
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Figure 14: Maximum eigenvalue of GJ^^^ versus the index distance d usin^ 
Regge action with A = cr = 1.0 (upper plot) and A = cr = 0.25 (lower plot). 
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Figure 15: Correlations between 4-volume elements using the Regge action with 
A = cr = 1.0 (upper plot) and A = cr = 0.25 (lower plot). 
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Figure 16: Largest eigenvalue ol the edge length correlation matrix versus d for a 
4-torus with 4"^ X 16 vertices using the Regge action. Two different cosmological 
constants A = 0.1, 1.0 are shown together with one coupling in the well-defined 
phase, /3 = 0.1, and with one coupling near the critical value, /3 = 0.14,0.108. 




Figure 17: Volume correlations on a 4-torus with 4"^ X 16 vertices. Cosmological 
constants A = 0.1, 1.0 are combined with /3 = 0.1 and with /3 = 0.14,0.108. 
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Figure 18: Curvature correlations on a 4-torus with 4"^ X 16 lattice. Cosmological 
constants A = 0.1, 1.0 are combined with different gravitational couplings. 
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6 Lattice Gravity as a Spin System 



It has been stated that the phase structure of simphcial quantum gravity shows 
certain similarities to spin glass models [29]. In this spirit, the path integral may 
be transformed to a partition function of a spin system using certain restrictions 
on the edges of the lattice [44]. This approach is structurally and computationally 
much simpler than the Regge calculus or the dynamical triangulation method. 
It is regarded as a third lattice gravity approach and is referred to as Ising-link 
quantum gravity in the literature [45]. 

6.1 Two dimensions 

The model is easy to define. All the squared edge lengths are allowed to take 
only two values 

qi = l + eai , < e < 0.6, a, G ^2 , (89) 

similar to the Regge-Ponzano model [46]. The real parameter e is restricted to 
fulfill the Euclidean triangle inequalities for the qi^s so that all configurations 
are allowed (A^i is the total number of edges). This is quite different from conven- 
tional Regge calculus where many potential updates either violate the triangle 
inequalities or the manifold property. Furthermore, it may provide a natural 
measure giving all 2^^ configurations equal weight. To use qi = k{l + eo";) is no 
more general, as k can be absorbed into a redefinition of A. 

According to the Gauss-Bonnet theorem the Einstein action in two dimen- 
sions is a topological invariant equal to Att times the Euler characteristic of the 
surface. Therefore, for a manifold with fixed topology only the cosmological term 
remains. To rewrite the action 

/ = A^Ai (90) 

t 

in terms of o"; we consider a single triangle (see Eigure 20). Its (squared) area 




Eigure 20: Notation for our triangular lattice. Triangle t consists of edges with 
the quadratic lengths qi^ q2, qi. 
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can be expressed as 
A] = 



Mqi + 92 - qi) 



0(91 + 92 - qi) 



92 



= + ^{cTi + (72 + cri)e + ^{aia2 + CTiai + a2cri - '^)e^ . (91) 

Expanding \fA^ the series consists only ol terms up to since af = 1. This 
suggests the following ansatz 

At = Co(e) + ci(e)(cri + (72 + (J/) + C2(e)(cricr2 + (JiCT/ + (72(7/) + C3(e)(7i(72(7/ . (92) 

There are only four possible values for the area of a triangle. Computing these ar- 
eas and comparing with (92) results in the following equations for the coefficients 

Ci 



Co 
Cl 
C2 
C3 



^[2y3 + 3/(e) + 3(/(e)] 
^[2V3e + fie)-gie)] 
^[2^3 -/(.)- (/(.)] 
^[2y3e-3/(e) + 3(/(e)] 



(93) 



where /(e) 



1 — e)(3 + 5e) and (/(e) 



1 + e)(3 — 5e), cf. Figure 21. Hence, 



one must have e < 7 for the triangle areas to be real and positive definite. Using 



Ci 



1 

0.8 
0.6 
0.4 
0.2 



-0.2 



1 1 1 1 1 
















_____ ^2!^— 


C3 ■ ■ ■ ■ ■ . 

1 1 1 1 1 



0.1 0.2 0.3 0.4 0.5 

e 



0.6 



Figure 21: Coefficients q [i = 0, ... ,3) obtained by rewriting the triangle area 
in terms of (7; = ±1. 
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(89) the measure (44) can be replaced by 



li [ ^ J2 exp[-m^ln(l + eai)] , 



(94) 



and the exponential in terms of e is reduced to 



m^ln(l + ecT/) = iVimo(e) + ^mi(e)cr/ , 



(95) 



with mo = — |me^ + 0{e'^) and mi = m(e + |e^) + 0{e^). 
Inserting (92) and (94) into the partition function yields 




t 



= J exp{-Y^[{2Xci + mi)ai + 



+ -Ac2(cri + (72 + 0-3 + (74)(7/ + -Ac3((7i(72 + (73(74)(7/]} , (96) 



with J = exp( — AA^2Co — Ninio) and A^2 the total number of triangles. Thus, 
the path integral becomes the partition function of a system consisting of a 
spin (7; at each edge /, with an external magnetic field and with 2- and 3-spin 
nearest neighbor interactions. Assigning the spin to the corresponding edge of 
the originally triangular lattice and drawing the interactions as lines a Kagome 
lattice is obtained (see Figure 22). 

Removing the term linear in cr by a convenient choice of the measure 
(mi = — 2Aci) and neglecting 3-spin couplings we get the partition function of 
an Ising model on a Kagome lattice. To compute the critical coupling we define 
Q = — |c2A and transform the partition function Zk{Q) of the Kagome lattice 
to the partition function Zd{L) of the decorated honeycomb lattice by applying 
the star-triangle transformation [47]: 



For (7i = (72 = (73 = 1 and (7i = (72 = —(73 = 1 we get 2cosh(3iv) = Aexp(3(5) 
and 2coshiv = Aexp( — Q), respectively. Any other choice of cr^ coincides with 
one of these two cases. Hence, A and Q are determined from 




Aexp[(5((7l(72 + (72(73 + (T^CTi)] 



^ exp[iv(7((7i (72 (73)] . 

(7 = ±1 



(97) 



A^ 

e4(3 ^ 



: e4Q(e4Q + 3)2 
2cosh(2iv) - 1 . 
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Figure 22: The triangular lattice is drawn by dashed lines and the Kagome 
lattice by solid lines. The spins are represented by o and • characterizing a 
typical configuration discussed in the text. 



Carrying out this transformation at every vertex gives Zd[L) = A~^^ Zk{Q) 
where Nh is the number of vertices of the honeycomb lattice. Furthermore, we 
can express Z^ by the partition function Zh of the conventional honeycomb 
lattice via the decoration-iteration transformation [47]: 



. L . L . — . . 

(Tl (7 (72 (Jl (72 

^ exp[ivcr(cri + (72)] = I ex]){K (71(72) . (99) 

(7 = ±1 

Similarly, we have 

/ = 2e^^' 

e2^^' = cosh(2iv) , (100) 
and inserting (98) and (100) we obtain 

ZoiL) = [e2^(e4^ + 3)]i^-ZK(g) = (2e^^')^-Z^(A0 

e^Q = 2e2^^' - 1 . (101) 

The critical couphng of the honeycomb lattice is given by e^^^'- = 2 + ^3 and 
thus for the Kagome lattice 

Qc = -|(Ac2)c = 0.4643 . (102) 
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In the ferromagnetic regime Q > (X < 0) the system shows a 2"''-order phase 
transition, whereas the antiferromagnetic regime Q < (A > 0) is governed by 
frustration even at zero temperature. 

Switching on the symmetry-breaking 3-spin interaction we performed nu- 
merical simulations of the system with toroidal topology in the coupling plane 
{Q2 = — |Ac2,(53 = "l^cs). Figure 23 shows the spin expectation value in a 
certain range of coupling constants. For Q3 = the conventional Ising model 
on a Kagome lattice is recovered and for 7^ 0, (^2 < the 3-spin coupling 
removes the frustration leading to ordered phases. Depending on sign((53) the 
spins denoted in Figure 22 by open dots take the values cr = ±1 whereas the full 
dots take cr = =pl giving rise to an expectation value (cr) = i|- 

Now return to the application of the spin system to the gravitational system. 
To investigate different types of measures the external field term has to be taken 
into account. Further the couplings depend on the value of e according to (93). 
Figure 24 depicts (cr) versus e for the uniform and the scale invariant measure 
and for the measure leading to a cancellation of the term linear in cr. The upper 
curve can be read off directly from Figure 23 and represents only a small detail 
near the origin. The 3-spin coupling is in the entire range of e not strong enough 
to remove the frustration and therefore (cr) 0. Both other curves show a 
favored occurence of negative spins, corresponding via qi = 1 -\- ecr; to short edge 




Figure 23: Spin expectation value as a function of the 2- and 3-spin couplings 
Q2 and Qs, respectively. 
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Figure 24: Spin expectation value versus e for nii = —2ci (no external field), 
m = (uniform measure) and m = 1 (scale invariant measure). In all cases the 
cosmological constant is set to A = 1. 



lengths reflecting the expected tendency of the lattice to shrink [48]. 

An extension of our approach to more than two edge lengths can be easily 
performed and it is straightforward to include additional scalar fields or physical 
spin fields. A generalization of the trivial 2-dimensional Regge- Einstein action 
to higher dimensions is of course possible, but becomes very complicated. Notice 
that in two dimensions only the three edges of a triangle are coupled in the 
action. In three dimensions one is faced with terms of 6^^ order at least, and with 
additional contributions from the Regge action. Recently, the Ising-link model 
in three dimensions has been analyzed via mean field techniques as well as with 
the Monte Carlo method [45]. However, critical exponents and the behavior of 
the curvature at the transition point differ from that found for Regge theory with 
continuously varying edge lengths. 

6.2 Four dimensions 

The situation is even more complicated in four dimensions where one has to deal 
with 10 edges in each simplex. Nevertheless, numerical simulations can be made 
very efficient by implementing look-up tables and using the heat bath algorithm. 
A spin update consists of choosing a particular spin, calculating the change in 
the action if the spin takes on its other possible value, and accepting the new spin 
value with probability to the exponential of the negative change in the action. 

In the actual computations we replaced (89) hy qi = bi[l -\- eai) because a 
4-dimensional Regge skeleton with equilateral simplices cannot be embedded in 
fiat space. Therefore, according to the indexing scheme of Section 3.1, hi takes 
different values depending on the type of the edge /. In particular hi = 1,2,3,4 
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Figure 25: Spin expectation value as a function of the gravitational coupling 
parameter for Ising-link quantum gravity in four dimensions. 
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Figure 26: Average action for simulations with different couplings fi and 
e = 0.0875. 



for edges, face diagonals, body diagonals, and the hyperbody diagonal of a hyper- 
cube, respectively. First preliminary studies show an interesting phase structure. 
The range of the gravitational coupling fi is extended also to negative values, be- 
cause there seems to be no particular reason to stay in the Fuclidean sector. 

The spin expectation value of a 4-dimensional toroidal lattice with 8"* vertices 
is depicted in Figure 25. The parameter e was set to 0.0875, A = and a = 1. 
The system seems to undergo phase transitions at two different points, first in 
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the negative region at /3 —4.7 and second at a large positive coupling /3 23. 

In Figure 26 the expectation value of the action versus the gravitational cou- 
pling /3 is presented. The resemblance of the curve for /3 > with those obtained 
by continuously varying edge lengths is very encouraging. One has to determine 
the value of the critical point and the order of the phase transition in the limit 
A^o — > oo in more detail. If both theories belong to the same universality class, 
one could perform calculations much easier and faster than with conventional 
Regge calculus. But to abandon dwelling on speculations, this is a question to 
be answered in future work. 
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7 Conclusion 



The Regge calculus provides a direct route to systematic approximations of the 
quantum-gravity path-integral. Unfortunately, there are almost no analytical 
results or experimental data that could guide us and thus, numerical simulations 
have to be interpreted very cautiously. In this spirit let us try to summarize. 

Although the Regge- Einstein action Ir and the "compact" action Ic are 
unbounded the entropy of the system stabilizes expectation values leading to 
the existence of a well-defined phase with small average curvature. Increasing 
the gravitational coupling /3 the system undergoes a phase transition, which is 
abrupt for Ir and smoother for Ic- This might be an important result for the 
continuum limit of simplicial quantum gravity. 

To get an idea of the details of the interaction, correlation functions between 
several observables like edge lengths, local volume elements, and local curva- 
tures have been computed. Most correlation functions exhibit a decrease of the 
corresponding effective mass towards the critical coupling. The infiuence of the 
cosmological constant shows no universal trend for the masses. However, no 
clear evidence for long range correlations is found, which should have indicated 
the occurence of massless particles. At present it is too early to preclude the 
existence of a massless graviton because with the appropriate action, measure, 
the true geodesic distance, and better statistics the reliability of the two-point 
functions will improve. 

Finally, a new method to investigate the phase structure of simplicial quan- 
tum gravity was introduced. Restricting the squared edge lengths of the lattice to 
take only two values the path integral can be transformed to the partition func- 
tion of a spin system. This allows to handle the system in a very economic way 
from an analytical as well as from a numerical point of view. Especially, if both 
quantum gravity and the corresponding spin system lie in the same universality 
class, this seems to be very promising for future investigations. 

If no 2"''^ order phase transition will be found, an intriguing idea would be to 
consider the lattice spacing as a fundamental length. Then one does not need 
to send the lattice spacing to zero [q — )• 0) as in usual lattice gauge theory, 
but there exists a limit qo = finite = 0{Lp). A further appealing idea is the 
construction of a grand unified theory of the fundamental forces by introducing 
additional matter fields, e.g. SU{2) gauge fields [49]. Then one has to try to 
tune the hadronic masses such that they are sufficiently small compared to the 
Planck mass. Also, the curvature of the lattice has to approach zero and the 
geometry to become fiat on the hadronic length scale. The infiuence of matter 
fields on the graviton mass would be of tremendous interest. 
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